Finite size and temperature effects in the J\ — J 2 model on a strip 
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Within Takahashi's spin-wave theory we study finite size and temperature effects near the quan- 
tum critical point in the Ji — J2 Heisenberg antiferromagnet defined on a strip (L x 00). In the 
continuum limit, the theory predicts universal finite size and temperature corrections and describes 
the dimensional crossover in magnetic properties from 2+1 to 1+1 space-time dimensions. 
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Takahashi's spin-wave theory is a useful tool for 
qualitative study of thermodynamic properties in low- 
dimensional spin systems u. For instance, the theory is 
capable to reproduce the universal low-temperature be- 
havior in two-dimensional (2D) Heisenberg antiferromag- 
nets (HAFM) □ originally predicted on the basis of the 
2+1 nonlinear a modcH in the limit p s <C T (p s is the 
spin- wave stiffness constant). Near the quantum critical 
point 2D HAFM is expected to possess universal mag- 
netic properties for arbitrary T/p s provided p s ,T <^ J, 
J being a characteristic exchange energy □. Finite size 
and temperature effects in 2D HAFM in the limit p s -CnT 
have previously been analyzed in the geometry L x LB. 

In the present work we apply Takahashi's approach to 
a frustrated Jj — J 2 Heisenberg antiferromagnet defined 
on a strip with linear dimensions L x 00 . In this geome- 
try the coupled self-consistent equations for the spectral 
gap A, the current field g — — (aQb'-), and the density 



field / = (a a±+y) can be represented in the for 
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where et = 4Jigr] k 1 y/l — T]%r/l is the magnon dispersion 
relation [rj^ 1 = 1 — (af/g){l — Kk)+fi, Kk = cosfci cosk 2 , 
7k = (cosfci +cosfc2)/2] a — J 2 /Ji is the frustration 
parameter, and S is the value of the lattice spin. The 
components of the wave vector k = fcix + k 2 y take the 
values fci = 2tttii/L [n x = 0, ±1, ... , ±(N - l)/2] and 
k 2 € [— 7r, 7r], where L is an odd integer. We use a system 
of units where h = kg = a = 1, a is the lattice spacing. 
The chemical potential p is connected to the spectral gap 
by the relation A = \ J\g\f2p . 

At T = and in an infinite geometry [L = 00) the sys- 
tem is in the magnetic Neel state, characterized by the 
on-site magnetization 



mo = S 
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up to the quantum critical point a c . Recent studies 
predict that the sysrem is magnetically disordered from 
a c = 0.38 up to a ps 0.6, and support the hypothesis that 
a c is a second-order phase transition point described by 
the 0(3) nonlinear a model 0. Eqs. ([jj) are supposed to 
describe the magnetically disordered phase at finite tem- 
peratures both in the infinite 00 x 00 and in the strip 
L x 00 geometries . In addition, the theory can also be 
applied for the magnetically disordered ground state of 
the strip system. In both cases the excitation gap A ^ 
appears as a result of the constraint that implies zero 
total staggered magnetization [the first of Eqs. (0)] and 
the theory is restricted to strip systems with the integer 
"rung" spin SL. An equivalent spin-wave description of 
the ground state of finite-size systems has originally been 
proposed by Hirsch and Tang. 

Leading corrections in T and l/L near the quantum 
critical point (T, a) = (0, a c ) can be obtained in the 
low-energy long-wavelength limit of the theory when the 
magnon dispersion relation takes the form 
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Here c = ^Jig(3 is the spin-wave velocity 

U3 = y/l/2 — otfjg^ and M = \/2p/ (3 is the magnon 

"mass" which is connected to the spectral gap by the 
relation A = cM. In the framework of the discussed 
theory the spin-wave stiffness constant is defined by the 
relation p s — 2Jig/3 2 mQ. In the long-wavelength limit 
Eqs. (1) are simplified to 
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In the last equations u) n = 2nnT and k = \Jk\ + k 2 . 
Finite-size effects at T=0: 
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M 2 <SW (L). (7) 



(5) For M « 1 the error function in the above equations 

can be replaced by unity and after some algebra one gets 
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and 



In the ground state the functions Wq{L) = W(L, 0) 
and Uo(L) = U(L,Q) take the forms 
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To extract the finite-size corrections we follow the 
method suggested in Ref. || (see also Re£jl0|). Using inte- 
gral representations for the integrandaliT/the above two 
equations are, respectively, transformed to 
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where erf(z) is the error function and Ql$) = 
(V^)E^ 2 : i 1 / / 2 2 + i /2 exp[-(27rn/L) 2 t]. Note that 

limi^oo Q L (i) = erf(7rv / t)/v / 47rt- 

The finite-size contributions in Wo(L) = Wo(oo) + 
SWo(L) and Uq{L) = Uq(oo) + SUq(L) can be separated 
by using-.the asymptotic formula (L » 1) (see Ref. [h], 
p. 148)0 
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The result reads 
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Here Li s (z) = J2T=i z "/ nS is the polylogarithmic function 
of s kind and 

F(y) = y 2 Lh(e- y ) + yU 2 ( c -y) + Li 3 ( e - y ). (10) 

Using Eq. (|^) , near the quantum critical point (T, a) — 
(0, q c ) the first of equations @) takes the form 
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Finally, a subtraction of the L = oo contributions in 
the last two of Eqs. (0) leads to the expressions 
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Here g c , / c , and f3 c are the values of the respective func- 
tions at the quantum critical point. The above expres- 
sions are valid near the quantum critical point. 

Eq. ( |Tl| ) exactly reproduces the saddle point equation 
(with an appropriate regularization) in the 1/N expan- 
sion of the O(N) nonlinear sigma model in 2+1 space- 
time dimensions. 113 As can be expected from the rel- 
ativistic dispersion relation (^), the ratio c/L plays the 
role of an effective temperature. It is well known that 
the gap equation (11) describes three different regimes, 
i.e., (i) the renormalizcd classical, (ii) the quantum crit- 
ical, and (iii) the quantum disordered regimes (see, e.g., 
Ref 0). 

(i) For p s L/c 1 (large linear-size behavior) the approx- 
imate solution reads A = Me = {c/L) exp(— 2np s L/(^. 
Using the mentioned symmetry and the refined resulto 
for the correlation length £ at low temperatures, the ex- 
pression for A can be improved to 
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where A is some dimensionless parameter. The last equa- 
tion may be interpreted as the spectral gap of a ID 
HAFM defined on a ladder-characterized by the integral 
"rung" spin S e ff = p s L/cM 

(ii) Exactly at the quantum critical point Eq. (O) gives 
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where yo = 21n(f/2 + y/E/2) is an universal constant. 
Eqs. ( |12] ) and ( |f3| ) lead to the universal finite-size 
contribution for the internal energy per spin [e(L) = 
-2J X ( 5 2 - a/ 2 )] 
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ttL 3 



F(LM), 



(16) 



where F(y) = F(y) — (y 2 /2)Lii(e _!/ ). It can be shown 
that for a < a c and in the limit LAI <C I the function 
F(LM) reduces to C(3), whereas exactly at the critical 
point we have F — 4((3)/5 + Uq/12: the last result is ob- 
tained by using Sachdev's identity (Ref. |l6|). Here £(x) 
is the Riemann zeta function. 

(iii)Finally, in the quantum disordered regime (a > a c ) 
we get 



e(L) = e c (oo) 
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The exponential form of the finite-size corrections reflects 
the existence of an excitation gap in the quantum disor- 
dered phase. Notice, however, that in this regime the 
system is known to possess gapped triplet excitations a, 
whereas Takahashi's theory is entirely based on trans- 
verse spin- wave excitations. 

Finite-size effects at low temperatures: 
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and the Poisson summation formula, one can transform 
Eqs. (§) into the forms W(L, T) = W(oo, 0) + 8W(L, T) 
and U[L,T) = f7(oo,0) + SU(L,T). In what follows we 
will consider only finite size and temperature corrections 
to the internal energy. Since at low temperatures the 
magnon mass M is exponentially small, it is enough to 
calculate the function SU(L,T) for M — 0. The result 
reads 
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\Ki(x) is the McDonald function] which yields the fol- 
lowing expression for the internal energy per spin: 
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In the considered long-wavelength approximation the ex- 
pansion parameter (TL/c)" 1 gives exponentially small 
contributions in the T 3 prefactor. In the bulk limit the 
T 3 term in Eq. ( po|) exactly reproduces the expected con- 
tribution to the internal energy of the 2D HAFM from 
the Goldg,tone modes at k = (0, 0) and k = (it, ir) (see, 
e.g., Ref.0). 

(ii) In the regime T -C c/L it is convenient to use in- 
stead of Eq. ( |19| ) the expression 
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which leads to the formula 
2<(3) 
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where the term m = n = must be omitted. 

(i) In the regime T c/L it is convenient to transform 

Eq. (|l8|) to the form 



L 3 ttcL 

In conclusion, based on Takahashi's spin-wave theory 
we have obtained expressions for the spectral gap and 
the internal energy per spin of a frustrated Heisenberg 
antiferromagnet defined on a strip. The expressions are 
valid in the continuum limit close to the quantum critical 
point and have a universal form depending only on the 
macroscopic magnetic parameters p s and c. 
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